Physical Chemistry Chemical Physics

Neutron scattering is a powerful technique to investigate atomic structures and molecular dynamics of proteins at the nano-scale. When it comes to dynamics, incoherent and coherent scattering respectively provide information on the single and collective dynamics of nuclei. In proteins, hydrogen has the highest incoherent cross-section, and it is common practice to overlook the contribution of coherent terms stemming from all nuclei. However, the fast collective dynamics of heavier nuclei could also be studied if coherent scattering and incoherent scattering were experimentally separated. The recent advent of polarized neutron spectroscopy with sufficient flux and energy resolution has made it possible, and opens new perspectives to investigate the relative importance of coherent scattering and the information it provides on biological samples. The present study reports on the use of polarized quasi-elastic neutron scattering (QENS) and the application of a minimalistic model adapted to both individual and collective dynamics. Using a perdeuterated green fluorescent protein as a model globular protein, the study provides an interpretation of the dynamical parameters obtained with QENS, and a comparative study of the Elastic Coherent and Incoherent Scattering Factor. Based on both experiments and calculations, we discuss the relative importance of distinct and self components of coherent scattering, which is often wrongly assumed to be representative of collective dynamics only. The results highlight the current impediments rendering complicated a straightforward analysis of fast collective dynamics in hydrated protein samples.

Adhesion of cells is the key factor determining functioning of multicellular organisms. Viscoelastic properties of cells can be studied by multiple methods. However, attractiveness of cells or extracellular matrix without the elastic component (dispersive adhesion) is not accessible. We present an extension of force spectrometry technology: the Multivalent Adhesive Probe Atomic Force Microscopy (MAPA) that delivers dispersive adhesion maps of live cells and biosurfaces, and identifies differences unresolved by viscoelastic probing.

We present the solution structures of HIV-1 proteins NC(p7)1-55 corresponding to the full-length NC(p7) and mature p6. The studies were carried in water and, to mimic the membrane, in micellar DPC (Dodecylphosphocholine) conditions. Our results unravel for the first time the structure adopted by the N-terminal amino acids of the free NC(p7)1-55, with the formation of a small helix spanning residues F6 to R10. Our NMR and Fluorescence Anisotropy data disclose an interaction between NC(p7)1-55 and p6 both in water and DPC, with respective Kd of 2.5mM and 370 mM at 23{degrees}C. The interaction is thus strengthened in lipidic conditions. Protein p6 stabilizes the N-terminus of NC(p7)1-55 while increasing at the same time the dynamic of the first zinc finger. Although the entire p6 sequence is involved in the interaction, we show that its C-terminal region is particularly sensitive to the presence of NC(p7)1-55, with a propensity of forming a a helix ranging from amino acids S111 to F116. This study brings experimental evidence of a direct protein-protein interaction between p6 and the N-terminal region of NC(p7)1-55. We further show that such interaction is readily accommodated within the NC(p15) framework and hypothesize that it may facilitate the selective assembly of assembly of the viral genomic RNA (gRNA) in the cell.

Carbapenems are the most potent {beta}-lactams, key antibiotics for healthcare-associated infections by Gram-negative bacteria and evade hydrolysis by most {beta}-lactamases, but are increasingly threatened by emergence of enzymes exhibiting hydrolytic activity towards them. Of the four recognised {beta}-lactamase subclasses, class A (active-site serine enzymes that hydrolyse {beta}-lactams via a covalent acylenzyme intermediate) is the most widely disseminated and, while the majority of such enzymes react with carbapenems to form long-lasting acylenzyme complexes, several possess carbapenem-hydrolyzing activity (carbapenemases). Here, we investigate the basis for these differences in a panel of class A {beta}-lactamases using molecular dynamics (MD) simulations of the respective acylenzyme complexes and tetrahedral intermediates (TI). The simulations reveal multiple features associated with catalytic activity across the spectrum of enzymes tested, including more extensive interactions of the carbapenem acylenzyme carbonyl and generally increased lifetimes of active site water molecules positioned for deacylation. Analysis of the dynamic trajectories shows carbapenemases to have reduced root mean-squared fluctuation (RMSF) differences between the acylenzyme and TI, that are not limited to the active site, indicating that the acylenzyme complex is pre-organised for reaction in carbapenemases but not in carbapenem-inhibited enzymes. Similarly, Principal Component Analysis (PCA) of acylenzyme and TI dynamics shows greater overlap between the two states in carbapenemases, providing further evidence for acylenzyme pre-organisation. Such simulations may represent an effective computational assay able to identify enzymes with carbapenemase activity at relatively modest computational cost.

Despite the well-known health risks of neutron exposures, key gaps remain in understanding neutron-induced molecular responses and identifying reliable biodosimetric markers that distinguish neutrons from photon exposure. We provide the first genome-wide analysis of the human blood transcriptional response to an accelerator-derived fission-like spectrum of neutrons versus photons, evaluating transcriptomic relative biological effectiveness (RBE) and radiation quality-discriminating gene signatures. Whole blood from healthy donors was irradiated ex vivo with X-rays (140 kV, 0-4 Gy, n = 3) or neutrons (0.1-8 MeV, 0-1 Gy, n = 2), incubated for 6 h or 24 h, and processed for RNA sequencing from peripheral blood mononuclear cells (PBMCs). Neutrons were markedly more potent than X-rays at inducing differentially expressed genes (DEGs) at equal doses, showing a peak response 6 h post-irradiation followed by a decline. In contrast, X-rays caused a continuous increase in DEGs up to 24 h (neutrons vs. X-rays at 1 Gy: 1,449 vs. 121 DEGs at 6 h; 996 vs. 621 DEGs at 24 h). A universal p53-centered 34-gene signature, including FDXR, EDA2R, GADD45A, and ZMAT3, showed highly monotonic dose responses (Spearman correlation coefficient {approx} 1) across donors, radiation qualities, and timepoints. Additionally, difference-in-differences analysis identified radiation quality-discriminating genes only at 6 h, with transcriptional convergence observed by 24 h, suggesting a very narrow time window for biodosimetric differentiation. We identified a neutron-specific gene signature driven by cGAS-STING-NF-{kappa}B signaling (RELB, NFKB1, C3, MALAT1) and suppression of B-cell and myeloid identity genes (IGHD, TCL1A, CLEC7A, TLR2), defining a biologically coherent neutron quality index with distinct immunomodulatory effects. For the first time, we assessed neutron RBEs at the gene, pathway, and global transcriptomic levels in a human blood model, reporting a global transcriptomic neutron RBE of 1.30 (95% CI: 1.14-1.49) at 6 h and 1.21 (95% CI: 1.14-1.28) at 24 h, providing a valuable basis for biodosimetry in mixed-field exposure scenarios. Our findings advance the mechanistic understanding of neutron radiation responses and support the development of biodosimetric approaches for mixed-field exposure scenarios.

Reading both physical and chemical properties of a microenvironment from a single fluorophore remains a challenge. Here we demonstrate that two coexisting molecular equilibria within one near-infrared cyanine, CyC4, encode two mechanistically distinct ratiometric reporting channels. A meso-amino group and a pendant carboxylate form a tunable intramolecular hydrogen bond that toggles the dye between closed (700 nm) and open (780 nm) emissive conformers. Time-dependent density functional theory (TD-DFT) calculations show that the hydrogen bond raises the LUMO and blue-shifts the emission, establishing the 700/780 emission ratio as a local reporter of hydrogen bonding and polarity. Independently, the chromophore self-associates under crowding- and cosolvent-rich conditions into an aggregate with a blue-shifted, H-type absorption signature near 530-540 nm and a distinct emission near 610 nm upon 540 nm excitation. The intensity of this aggregate band relative to the monomer emission (Ra) serves as a ratiometric reporter of crowding and self-association. Because the two channels arise from distinct molecular equilibria (intramolecular hydrogen bonding vs. intermolecular self-association) they are largely decoupled: a glycerol titration series confirms that the self-association channel (Ra) can be moved while the hydrogen-bonding channel stays essentially fixed. Applied to protein-PEG biomolecular condensates, the two ratios move oppositely with increasing salt, showing that the interior's chemical (polarity, hydrogen bonding) and physical (packing, self-association) environments co-vary across the salt series; a single CyC4 measurement thereby maps this coupled microenvironment, providing a general strategy for multiparametric, ratiometric sensing of crowded microenvironments.

This paper presents and critically reviews the results of original and some literature based experimental studies conducted by the authors last years on the structural organization of DNA in dormant (starvation stress), anabiotic dormant (4 HR treatment) E. coli cells, as well as the K12 {Delta}dps strain, which lacks the Dps protein (Dps null E. coli). The experimental data includes small-angle synchrotron radiation diffraction (SAXS) and transmission electron microscopy (TEM) data. Synchrotron radiation diffraction experiments on K12{Delta}dps cells allowed us to conclude that peaks at 44.3, 22.1, and 14.8 angstrom resolutions are associated exclusively with ordered DNA organization. Peaks at 44.3, 22.1, and 14.8 angstrom resolutions are also observed for samples of dormant (starvation stress) cells and anabiotically dormant cells. Therefore, this ordered DNA organization also applies to samples of dormant and anabiotically dormant cells. A model is proposed that considers the ordered DNA organization in the cell as a cholesteric liquid crystal. The powder diffraction pattern calculated based on this model is compared with experimental small angle X ray scattering (SAXS) data obtained on Dps-null cell samples. The model completely reproduces the key features of the experimental diffraction pattern from Dps-null cell samples. Accordingly, the cholesteric liquid crystal model corresponds to DNA packaging in dormant and anabiotically dormant cells. Cholesteric liquid crystal ordering should be further considered in all models of cellular DNA packaging. To address the question of which structural organization of DNA predominates in the cell: the cholesteric liquid crystal or nanocrystalline or whether they coexist and fully manifest themselves under different external conditions, it is necessary to utilize the latest methodological advances in structural analysis.

UCH-L1 is a monomeric deubiquitinating enzyme whose native structure embeds a shallow $5_2$ knot located near the N-terminus, placing the knotted topology in direct proximity to both the substrate-binding pocket and the catalytic site. While our previous work established that N-terminal integrity is critical for catalytic activity, the energetic cost of unknotting and its structural consequences remained unquantified. Here, we combine steered molecular dynamics with an umbrella sampling scheme to generate topologically modified variants of UCH-L1 and, for the first time, reconstruct the free-energy profile of UCH-L1 unknotting. The potential of mean force reveals a steep energetic barrier to knot disruption, consistent with knotting being a late, rate-limiting folding step that is effectively locked in once the native structure is established. Long unbiased MD simulations of fully unknotted variants in both apo and holo states show that knot removal increases local flexibility at the N-terminus without inducing significant global structural destabilization. Binding energy calculations indicate that the unknotted variant binds to ubiquitin less tightly than the wild-type ($\sim$-62~vs~$\sim$-76~kcal/mol), suggesting that topological integrity contributes to substrate affinity. Together, these results show that the $5_2$ knot in UCH-L1 is not a passive structural feature but a functional element that fine-tunes folding kinetics and contributes to substrate binding efficiency.

Current approaches for imaging fibrotic remodeling have sensitivity, specificity and cost drawbacks that limit both preclinical research and clinical diagnosis. Here, we show that fast green FCF, a small molecule that binds to fibrillar collagen, enables highly sensitive and specific imaging of fibrosis in lung samples from mice and humans using fluorescence microscopy. We report strategies for using fast green FCF staining to assess fibrotic remodeling using precision-cut lung slice and whole-biopsy preparations. Our findings demonstrate that fluorescence imaging of fast green FCF-stained collagen will be useful for fibrosis research and may help to improve detection of fibrosis in clinical pathology.

Online structural characterization during microfluidic lipid self-assembly is important for understanding and controlling the formation of nonlamellar liquid crystalline nanodispersions. Here, we report a 3D-printed, X-ray-compatible hydrodynamic flow-focusing microfluidic chip with variable channel dimensions, integrated with synchrotron small-angle X-ray scattering (SAXS), for position-resolved SAXS-on-chip monitoring of Ca2+-triggered hexosome formation. Hexosomes were produced under continuous flow by mixing ethanolic solutions of docosahexaenoic acid monoglyceride (MAG-DHA), the negatively charged phosphatidylglycerol DOPG, and -tocopherol with Ca2+-containing PIPES buffer. Online SAXS-on-chip measurements detected three Bragg reflections characteristic of the internal inverse hexagonal (H2) phase on a tens-of-milliseconds residence-time scale, revealing rapid structural evolution during microfluidic mixing. Complementary ex situ SAXS identified the DOPG/Ca2+ molar ratio as a key parameter modulating the direct vesicle-to-hexosome transformation and the compactness of the internal H2 nanostructures. Dynamic light scattering showed that the flow-rate ratio modulated nanoparticle size, yielding hexosomes with mean hydrodynamic diameters in the range of approximately 120-175 nm and polydispersity index values down to 0.14 at a total flow rate of 200 {micro}L min-1. Cryo-TEM revealed coexistence of hexosomes and vesicular nanostructures, highlighting morphological heterogeneity, while Coarse-Grained Molecular Dynamics simulations supported a central role of Ca2+-DOPG association in promoting a direct lamellar-H2 phase transition. Overall, this work shows that 3D-printed SAXS-compatible microfluidics can integrate continuous production with online structural characterization, providing a basis for future formulation and process optimization of drug-loaded cubosomes, hexosomes, and related nonlamellar liquid crystalline nanodispersions.

The navigational prowess of migratory birds is thought to arise from light-dependent radical-pair chemistry in cryptochrome 4 (CRY4), yet the slow structural transitions that couple photochemistry to signaling remain elusive. Here, we combine temperature-controlled steady-state UV-visible spectroscopy and hydrogen-deuterium exchange mass spectrometry (HDX-MS) to elucidate the photochemical and conformational dynamics of pigeon CRY4 (ClCRY4). Steady-state measurements at 5-25 {degrees}C reveal that lower temperatures slow FAD photoreduction and prolong the FAD neutral semiquinone signaling state. This occurs without a solvent kinetic isotope effect, implicating a conformational change rather than proton transfer as the rate determining step in FAD neutral semiquinone formation. Simultaneous HDX-MS under blue-light exposure identifies protection near the FAD-binding site and C-terminal region. To enhance sensitivity, we developed a pump-probe HDX-MS approach at 10 {degrees}C. This reveals eight peptides (within the phosphate-binding loop, protrusion motif, electron-transfer-chain loops and C-terminal tail) that exhibit rapid ([≤]10 s) and sustained light-induced protection, delineating early conformational rearrangements as a prerequisite for FAD neutral semiquinone accumulation. The findings of slower onset HDX protection as well as a bimodal pattern of deuterium uptake in the phosphate-binding loop further identify a local redistribution of conformational substates on the time scale of the accumulation of the signaling species. Site specific mutagenesis within the CTT supports the findings, which lead to a model in which blue light triggers rapid clamping down of protein near the two regions of spin pair separation, followed by a rate limiting closure of a surface loop. The resolution of time-dependent structural transitions that follow photoactivation of CRY4 resolves the interface between quantum radical-pair formation and classical conformational changes, while providing an enhanced structural framework for the molecular events that underlie avian magnetoreception.

Cyclic oligosaccharides possess industrial advantages, including molecular encapsulation capability and high physicochemical stability, owing to the absence of a reducing end. Recently, a novel cyclic tetrasaccharide, cycloisomaltotetraose (CI4), consisting of four -1,6-linked glucose units, and the enzymes responsible for its synthesis, cycloisomaltotetraose glucanotransferases (CI4Tases), were discovered. Unlike known cycloisomaltooligosaccharide glucanotransferases (CITases) that yield a wide distribution of cyclic products with a degree of polymerization (DP) of 7 or higher, CI4Tases strictly produce CI4. To elucidate the molecular mechanism underlying this strict DP4 specificity, we determined the crystal structures of CI4Tase from Agreia sp. D1110, in its ligand-free form, as well as in complex with the linear hydrolysis product isomaltotetraose (IG4) and with CI4. Structural comparisons revealed that a loop (M247 to R251) blocks the region corresponding to the -5 subsite of typical CITases, narrowing the substrate-binding pocket. This "molecular ruler" mechanism ensures that only a glycan chain of exactly four glucose units is accommodated for cyclization. Among mutants of the residue positioned at the center of bound CI4, the formation of by-products other than CI4 was significantly suppressed in F245L, F245A, and F245W. While the cyclization activity of all F245 mutants decreased, the CI4 hydrolysis activity of these three mutants was also significantly reduced, resulting in an increased specificity for cyclic sugar production. These findings elucidate the strict size-control mechanism of CI4Tase and provide a structural foundation for engineering cycloisomaltooligosaccharide-producing enzymes with optimized transglycosylation efficiency and specificity for industrial applications.

In linear polysomes, excluded-volume interactions among ribosomes can induce dimensional reduction of mRNA. Yet linear architectures allow steric stress to relax at open ends-- limiting how strongly crowding can remodel the mRNA's structure and dynamics. Using coarse-grained molecular-dynamics simulations, we compare circular and linear polysomes over a range of ribosome densities. Circular closure selects a predominantly quasi-planar global conformational ensemble, as indicated by a shape dimensionality dshape {approx} 2 over a range of ribosome densities. Crucially, circular topology and ribosome crowding act cooperatively to suppress structural fluctuations. While closure alone or linear crowding reduces relative global size fluctuations ({Delta}Rg/Rg) only to {approx} 0.16, their combined effect drives this fluctuation down to {approx} 0.07. Within this stabilized architecture, increasing ribosome density drives a distinct in-plane reorganization: the ring becomes more isotropic, global size fluctuations are strongly suppressed, and the scaling exponent increases toward {nu} [~=] 0.74 - 0.77, consistent with two-dimensional self-avoiding walk-like value over the accessible finite-size window, 1000 [≤] N [≤] 4969. Closure shortens the radius-of-gyration decorrelation time of circular polysomes by 40-fold relative to matched linear systems, reflecting the topological elimination of free ends. Within this closureselected ensemble, ribosome crowding further reduces the decorrelation time by up to 20% at the highest density. A fluctuation-informed crossover model links the density dependence of the global scaling exponent to inter-ribosomal subchain statistics. These results distinguish the geometric role of circular closure from the density-dependent steric response that it enables, revealing a confined yet dynamically responsive conformational regime for circular polysomes.

Many biological processes rely on mechanical forces, with protein molecules acting as key mediators. Understanding how proteins respond to mechanical stress is essential for conditions including cardiomyopathy and muscular dystrophy. Natural proteins such as dystrophin and utrophin are composed of heterogeneous folding domains with distinct mechanical properties; deciphering domain-level behavior provides insights into disease mechanisms and informs therapeutic strategies. Single-molecule force spectroscopy (SMFS) enables probing the mechanical properties of entire proteins, yet current approaches struggle to identify heterogeneous folding domains, particularly without prior knowledge. Here, we present the first automated framework to identify heterogeneous folding domains in SMFS data, applying both existing clustering methods and a novel physics-aware deep clustering architecture, LatentUnfold. LatentUnfold learns complementary latent representations from force magnitude and the force-extension physical relationship through dual autoencoders, jointly optimized for clustering assignments. We apply our framework to experimental SMFS data collected from a synthetic two-domain protein (ddFLN4-Titin I27) as well as natural protein constructs of dystrophin and utrophin, with Monte Carlo simulated datasets serving as controlled validation. For the synthetic protein, we recover mechanical properties consistent with previously reported values for each domain. For the natural proteins, we uncover two mechanically distinct domain populations - corresponding to the N-terminal domain and spectrin-like repeats - with differences in both unfolding force and contour length increase, and reveal different unfolding order between them for the first time. This work enables domain-level biological inference, overcoming prior limitations that relied on averaging and overlooked heterogeneity, thus advancing the understanding of mechanical behavior in protein unfolding.

The binding affinities of a chimeric analog of a riboflavin derivative linked to biotin, (6- (7,8-dimethyl-2,4-dioxo-3,4-dihydrobenzo[g]pteridin-10(2H)-yl)hexyl 5-((3aS,4S,6aR)-2- oxohexahydro-1H-thieno[3,4-d]imidazol-4-yl)pentanoate), referred to as C6-Rf-biotin-tag, to the riboflavin binding retain or streptavidin are in the M range, 1.29 {+/-} 0.277 and 3.00 {+/-} 0.459, respectively. These values suggest that C6-Rf-biotin-tag has potential applications in diagnostic assay and labelling target flavin binding proteins. The C6-Rf-biotin-tag which was characterized with respect to physical and biochemical properties retains UV/Vis spectroscopic and fluorescence behavior similar to riboflavin.

Transcription factor Ets2 coalesce with the NF-{kappa}B pathway to regulate gene expression in specific signaling contexts. IKK2/{beta}-mediated phosphorylation events critically regulate the NF-{kappa}B pathway. However, any link between Ets2 and IKK2 remains elusive. Here we report Ets2 as a direct substrate of IKK2. In-vitro kinase assays using deletion constructs, high resolution MS-MS and site directed mutagenesis identified S295 as a prominent phosphorylation site distal to the DNA binding domain, substitution of which to phosphor-mimetic Glutamate triggers further phosphorylation of Ets2. MD simulations clearly indicate conformational constriction of the otherwise disordered N-terminal region and inhibition of DNA binding activity upon phosphorylation, which was further confirmed by Electrophoretic mobility shift assays. Our results uncover a phosphoregulatory connection between Ets2 and IKK2.

Amyloid fibrils are highly ordered protein aggregates characterized by a conserved cross-{beta}-sheet architecture despite originating from structurally diverse precursor proteins. Growing evidence suggests that interactions between different amyloidogenic proteins can modulate aggregation pathways through heterologous cross-seeding; however, the influence of seed polymorphism on the structure and biological properties of cross-seeded fibrils remains poorly understood. Here, we investigated the cross-seeding of native human insulin by two structurally distinct polymorphs of hen egg-white lysozyme (HEWL): flexible fibrils (FFs) and rigid fibrils (RFs). Native insulin remained stable under physiological conditions and underwent spontaneous fibrillation only under acidic conditions. In contrast, both HEWL polymorphs efficiently induced insulin aggregation at physiological pH, bypassing the nucleation barrier. Thioflavin T fluorescence, circular dichroism spectroscopy, and transmission electron microscopy revealed that lysozyme FFs templated the formation of insulin flexible fibrils (IFFs), whereas lysozyme RFs produced insulin rigid fibrils (IRFs), demonstrating that the structural characteristics of the parental HEWL polymorphs were propagated during heterologous cross-seeding. The toxicity of the resulting insulin fibrils was evaluated in SH-SY5Y neuronal cells and CCF-STTG1 astrocytes. IFFs exhibited minimal cytotoxicity and only subtle morphological alterations, whereas IRFs caused modest reductions in cell viability accompanied by more pronounced cellular damage. These findings demonstrate that the structural polymorphism of HEWL fibrils governs both the architecture and biological activity of cross-seeded insulin fibrils, highlighting amyloid polymorphism as an important determinant of heterologous amyloid propagation and a potential design principle for engineering functional amyloid-based biomaterials and protein delivery platforms.

Protein synthesis supports most biological processes. In the brain in particular, protein synthesis plays a critical role in physiological and pathological states. Here, we describe Tellurophene-Alkyne Cycloaddition-mediated Amino acid Tagging (TeACAT), a versatile strategy for fast, facile, and flexible tagging of newly synthesized proteins in mice. TeACAT is based on metabolic incorporation of the non-canonical amino acid TePhe into proteins by the endogenous protein synthesis machinery. Due to their high similarity, TePhe can efficiently replace canonical Phe without dietary or genetic manipulation. The subsequent bio-orthogonal reaction of TePhe with either fluorescent dyes or affinity handles enables both visualization and affinity enrichment of proteins synthesized during TePhe exposure. TeACAT is compatible with immunofluorescence for cell-type specific visualization of protein synthesis with subcellular resolution and can be used in conjunction with routine proteomics to identify and quantify newly synthesized proteins. Robust incorporation into the mouse proteome was observed on the scale of hours to days, allowing the interrogation of various biological processes. In summary, TeACAT enables the visualization and quantification of protein synthesis with minimal perturbation for biological discoveries.

The folding thermodynamics of proteins are dominated by two opposing forces, the loss in backbone entropy and the packing of hydrophobic groups. The same forces are major contributors to the extension thermodynamics of elastic proteins with the distinction that both processes act in concert, favoring the higher chain and solvent entropy of a relaxed conformation. The relative entropic contributions specify the recoil mechanism; human elastin recoil is primarily driven by hydrophobic forces, whereas fly resilin has a rubber-like mechanism driven by backbone entropy. Despite the importance of elastic proteins to tissue biomechanics, few have been identified, let alone characterized to the same extent as elastin and resilin. We develop a thermodynamic framework that maps proteins by sequence-derived estimates of extension-induced backbone and solvent entropy changes. Putative elastic proteins are proposed and classified by recoil mechanism based on estimated thermodynamic features. Proteins that map to elastic regions are overrepresented by the skin proteome. The set of predicted elastic domains is further extended by incorporating sequence context embedded in protein language models. Protein domains with distinct thermodynamic recoil mechanisms cluster on the latent space manifold. Some of these domains are anticipated to have roles within molecular machines, expanding the scope of elastic protein function beyond mechanical materials like elastin and resilin.

The cluster of differentiation 36 (CD36) is a membrane protein with broad physiological roles in health and disease, and its function is regulated in part by phosphorylation. Experimental evidence shows that phosphorylation of Thr92 reduces CD36 affinity for thrombospondin-1 (TSP-1), binding of which initiates antiangiogenic signaling, whereas phosphorylation of Ser237 decreases CD36-mediated fatty acid uptake, with implications for energy metabolism. However, the only available crystal structure of CD36 lacks phosphorylation, and the molecular mechanisms by which phosphorylation regulates CD36 function remain largely unknown. This study provides an atomically detailed computational characterization of CD36 in unphosphorylated and dual phosphorylated states, using molecular dynamics simulations with a total sampling time of 30 microseconds in combination with Markov state models. We present, to our knowledge, the first evidence of a cryptic pocket on CD36 surface that is formed by phosphorylation. This cryptic surface pocket and a loop spanning residues 121-131 form a high affinity binding site for TSP-1 derived ligands, shifting their binding away from the canonical site. We propose that this altered binding provides a molecular basis for the disruption of antiangiogenic signaling upon CD36 phosphorylation. Additionally, our data indicate that, phosphorylation increases helicity and compaction within the helix-loop region spanning residues 296-331, narrowing one of the entrances to the internal cavity and reducing its overall volume. These conformational changes provide a potential mechanistic explanation for the decrease in fatty acid uptake upon CD36 phosphorylation. Our findings provide structural insights that may inform the future design of CD36 modulators and emphasize the importance of targeting phosphorylation induced CD36 conformations in angiogenic and metabolic diseases.